NTOUW-02-002 
Nucleon-Nucleon Interactions on the Lattice 



Silas R. Beane and Martin J. Savage 

Department of Physics, University of Washington, 
Seattle, WA 98195. 



Abstract 

We consider the nucleon-nucleon potential in quenched and partially- 
quenched QCD. The leading one-meson exchange contribution to the po- 
tential is found to fall off exponentially at long-distances, in contrast with 
the Yukawa-type behaviour found in QCD. This unphysical component of the 
two-nucleon potential has important implications for the extraction of nuclear 
properties from lattice simulations. 



February 2002 



1 



It is hoped that lattice methods will ultimately lead to an understanding of nuclear 
physics directly from QCD. To date there exist two distinct approaches to investigating the 
nucleon-nucleon (NN) system using lattice QCD. The first approach makes use of Luscher's 
finite-volume algorithm, which expresses the energy of a two-particle state as a perturbative 
expansion in the scattering length divided by the size of the box |L|. Practical computations 
require a sufficiently small box, which in turn should be larger than the scattering length. 
Hence Luscher's method is ideal for systems with "natural" scattering lengths - with size of 
order a characteristic physical length scale - and has been used to study im scattering 0. 
Of course in the NN system, the scattering lengths are much larger than the characteristic 
physical length scale given by A^^. Nevertheless one may expect that at the unphysical 
values of the pion mass currently used in lattice simulations, the scattering lengths relax 
to natural values, thus allowing their determination from the lattice. Attempts have been 
made to compute NN scattering parameters in lattice QCD; Ref. || computes the : 5o and 
3 Si scattering lengths in quenched QCD (QQCD) using Luscher's method. A second, less- 
ambitious approach to the NN system on the lattice is to study the simplified problem of 
two interacting heavy-light particles |§-|8|]. Here the hope is to extract information about 
the long-distance tail of the potential, which in nature is provided by one-pion exchange 
(OPE). 

At present computational efforts are limited by the use of quark masses that are sig- 
nificantly larger than those of nature, typically rn}^ u - ~ 500 MeV. In order to make a 
connection between present (and near-future) lattice calculations and nature, an extrapola- 
tion to the physical values of the quarks masses is required. The technology which allows a 
systematic extrapolation has been developed for QQCD P-p~3|j and partially quenched QCD 
(PQQCD) []I4|-|I7|. In this paper we will consider the potential between two nucleons at lead- 
ing order (LO) in quenched chiral perturbation theory (Q%PT) and in partially-quenched 
chiral perturbation theory (PQ^PT), the low-energy effective field theories that describe 
QQCD and PQQCD, respectively. 

We will be considering the extension from two-flavor QCD, where the three 7r's are 
pseudo-Goldstone bosons associated with the breaking of chiral symmetry and the 77 is the 
SU(2) singlet meson whose mass arises mainly from the U{1)a axial anomaly, to QQCD and 
PQQCD. Essential to our discussion will be the modifications to the r\ propagator which 
result from quenching and partial quenching. In QQCD the quarks, Q = (u,d,u,d) T , live 
in the fundamental representation of SU (2\2), with the bosonic quarks, u, d, accounting for 
the graded structure of the symmetry group. The low-energy effective field theory (EFT) is 
QxPT, whose degrees of freedom are constrained by the chiral group SU(2\2) L Cg> SU(2\2) R . 
In the isospin limit the rj propagator in Q^PT is given by [0-|TTJ 

(q z — + is) {q z — m~. + ie) z 

The hairpin parameters M and a® characterize the dynamics of the singlet field and do not 
vanish in the chiral limit. The double-pole, or hairpin, structure is an unphysical artifact 
of quenching. Unfortunately QQCD is not related to QCD by any smooth variation of 
parameters. This is not the case in PQQCD where the quarks, (four fermionic and two 
bosonic) Q = (u, d,j, I, u, d) T live in the fundamental representation of the graded symmetry 
group SU (4|2) and the low-energy EFT is PQxPT, whose degrees of freedom are constrained 
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FIG. 1. Pion and r/ exchange contributions to the NN potential. The cross signifies a hairpin 
propagator from eq. ([|) (QQCD) or eq. (PQQCD). 



by chiral SU(4\2)l <8> SU(4\2)r. As the j, I sea quarks become degenerate with the u, d 
valence quarks, PQQCD smoothly maps to QCD. If the sea quark masses are taken to 
infinity, PQQCD formally goes to QQCD. In the isospin limit for both valence and sea 
quarks the r\ propagator in PQ^PT is 



where m Vj is the mass of the pseudo-Goldstone bosons composed of sea quarks. Clearly 
this unphysical double pole vanishes in the QCD limit. The hairpin structure resulting from 
(partial) quenching will have important implications for lattice simulations of heavy-heavy 
systems, such as the two-nucleon system or systems composed of two (or more) hadrons 
containing heavy quarks. 

The EFT power-counting for NN scattering has been developed in Ref. [|T§| -|20||. While 
details of the power-counting are channel dependent, LO for the low partial waves consists 
of OPE and a four-nucleon contact interaction with no derivatives or insertions of the quark 
mass matrix. This contact operator encodes information about short- distance physics. The 
EFT power-counting is readily generalized to QQCD and PQQCD, and there is no direct^ 
(partial) quenching of the NN potential at LO in the EFT. This is easily seen by assigning 
a conserved charge to the various quark types. We will consider PQQCD but the same 
argument applies to QQCD. For instance, consider the charge assignments qy = (1/3,0,0), 
qs = (0,1/3,0), qc = (0,0,1/3) for valence-, sea- and ghost-quark number, respectively. 
Here qy plays the role of ordinary baryon number. The NN initial state carries total quark- 
number charge (2, 0, 0). Baryon states with two valence quarks and either one ghost quark or 
one sea quark carry charge (2/3, 0, 1/3) or (2/3, 1/3, 0), respectively, while meson states with 
one valence quark and either one ghost or one sea anti-quark carry charge (1/3, 0, —1/3) or 
(1/3, —1/3, 0), respectively. It follows that the NN potential must have NN in the initial and 
in the final state. Of course the potential is directly (partially) quenched, but only through 
loops, which appear beyond LO in the EFT, as in fig. (^|b). Direct (partial) quenching 
modifies the tree-level potential through higher fock components. For instance, one can 



1 We define direct (partial) quenching to be (partial) quenching beyond modifications to the singlet 
propagator. 
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FIG. 2. The solid lines indicate valence quarks and the dashed lines indicate either a sea quark 
or a ghost quark. In the left panel (a) we exhibit a higher fock component of the NN potential and 
in the right panel (b) we exhibit a crossed-box loop correction to the NN potential. Both graphs are 
subleading in the EFT expansion. 

have a final state of one nucleon and one baryon composed of two valence quarks and one 
sea or ghost quark and one meson composed of one valence quark and one sea or ghost anti- 
quark. This contribution can arise either from meson-exchange graphs, as in fig. (0a), or from 
four-nucleon contact operators with one insertion of the quark mass matrix or one insertion 
of a chiral covariant derivative. In either case, the external meson legs at the level of the 
potential must be reabsorbed by the final state nucleons. In contrast with the "potential" 
pions exchanged between the baryons, these "radiation" pions generate contributions which 
are subleading in the EFT. Hence at LO there is no direct (partial) quenching and the 
NN potential is given by the one meson exchange graphs of fig. (0). However, there are 
hairpin contributions to the potential at LO. For simplicity we work in the isospin limit for 
both the valence and the sea quarks. The diagrams in fig. (|l|) give the LO meson-exchange 
contributions to the NN potential. In coordinate space the NN potential is 

V™(r) = 1 n . Va 2 ■ V U ^ + gS - si M °: at '"- ) e — , (3) 

oirj z \ r r 2m n J 

in QQCD, where the axial coupling to the pions is g^, while the coupling to the rj is go- 
However, one should keep in mind that there is no relation between axial couplings in QQCD 
and those of QCD. In PQQCD, we find that 

v™ M = ^ * ■ v* • v fa ^ - B l M^l) e ~ , ,4) 

where the couplings qa and go in this potential should be perturbatively close to those of 
QCD. Notice that effects due to partial-quenching (but not quenching) disappear in the 
short-distance limit, and are irrelevant to the renormalization of the potential. This in turn 
implies that the four-nucleon contact operator, that appears in PQQCD at the same order 
in the EFT as OPE, is that of QCD |T8^pO|. In contrast, the long-distance behaviour of the 
potentials are dominated by the hairpin contributions, 
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(m^ -ml) , PQQCD 



(5) 



which has both a central piece and a tensor piece. Both quenching and partial-quenching 
lead to an NN potential which falls off exponentially at long distances, as opposed to the 
physical Yukawa fall-off. In nature there is a fine-tuning between the long-distance Yukawa 
part of the potential due to OPE, and the short-distance part which is encoded in the contact 
operator which appears at the same order in the EFT. It is this fine-tuning which leads, for 
instance, to a shallowly bound deuteron. It is clear that this fine-tuning does not exist for 
QQCD, and quickly disappears in PQQCD as one moves away from the QCD limit. 

In conclusion, multi-nucleon quenched and partially-quenched simulations on the lattice 
are complicated by the presence of a long-distance exponential tail, that disturbs a fine- 
tuning which is a key feature of nuclear physics. It has been shown in Ref. [OT that an 
understanding of the quark-mass dependence of the deuteron binding energy requires the 
matrix element of an operator that is not constrained by NN scattering data but which in 
principle can be obtained from lattice QCD. Any meaningful determination of this matrix 
element from a partially-quenched simulation will have to account for the unphysical, long- 
distance behavior highlighted in this work. 
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